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Abstract 

In this paper is discussed an application of signed measures (charges) to description of 
segment and chord length distributions in nonconvex bodies. The signed distribution may 
naturally appears due to definition via derivatives of nonnegative autocorrelation function 
simply related with distances distribution between pairs of points in the body. In the work is 
suggested constructive geometrical interpretation of such derivatives and illustrated appear- 
ance of "positive" and "negative" elements similar with usual Hanh-Jordan decomposition in 
measure theory. The construction is also close related with applications of Dirac method of 
chords. 
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1 INTRODUCTION 



Different properties of chord length distribution (CLD) for nonconvex bodies were discussed re- 
cently in few publications [U [H [31 [H [31 [7] . Let us recall three different ways to introduce CLD 
for nonconvex body. Straight line may intersect nonconvex body more than one time (see Fig. [1]) 
and we can either consider each segment of such line as separate chord or calculate sum of all such 
segments. These two methods are known as multi-chord and one-chord distribution (MCD and 
OCD) respectively [Il[l[3]. 




Figure 1: Nonconvex bodies: a) simply connected b) with hole. 

For the convex case probability density function for CLD is proportional to second derivative 
of autocorrelation function [U HI |5l H] . Such a property may be used for a third definition of "the 
generalized chord distribution" [31 [5l [6] , but straightforward calculations demonstrates possibility 
of negativity of such function for some nonconvex bodies [6] . Let us call this function here signed 
chord (length) distribution to avoid some ambiguity of term "generalized" and to emphasize the 
basic distinguishing property of this function. 

Definition of CLD for convex body has standard probabilistic interpretation in theory of geo- 
metric probability and random sets [8l[9l[T0]. The MCD and OCD cases for nonconvex body may 
be described as well [2 [21 [3] . Is it possible to consider similar possibility for the signed chord 
distribution! 

In wonderful essay "Negative probability" [llj Feynman wrote that, unlike "final probability of 
a verifiable physical event" , "conditional probabilities and probabilities of imagined intermediate 
states may be negative" and so: "// a physical theory for calculating probabilities yields a negative 
probability for a given situation under certain assumed conditions, we need not conclude the theory 
is incorrect. " In this review Feynman provided a few examples with appearance and interpretation 
of negative probabilities both for quantum and classical physical models. 

Mathematical extension of the measure theory for such a purposes may use so-called signed 
measures ( charges ) |12j . Usually such extension is reduced to standard positive measures due 
to Hahn and Jordan decompositions, corresponding to expression of charge as difference of two 
positive measures [12] . 

In many processes with signed distributions the Hahn decomposition, i.e., splitting of space 
of events on positive and negative parts is quite obvious, e.g., in simplest examples we have two 
kinds of events: putting and removing objects [11]. A distinction of the signed chord distribution 
is appearance of negativity due to differentiations of positive function without such a natural 
decomposition on positive and negative elements. 

Nonconvex body Fig.[T|D may be represented as a convex body and a convex hole and it provides 
some intuitive justification of possibility to express some distributions using formal difference of 
convex hull and the hole. Rigor consideration is more difficult, especially for chord distribution 
expressed via second derivative, e.g., method derived below in Sec. l3.3l reduces examples like Fig.[T] 
to six "signed" intervals: four "positive": [AD], [AB], [CD], [BC] and two "negative": [AC], [BD]. 
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The convex case is revisited in Sec. [5] and Appendix |^ The construction of signed chord 
length distribution for nonconvex body is described in Sec. [3] Some imphcations to description 
of arbitrary bodies with nonuniform density are briefly mentioned for completeness in Sec. 2] and 
Appendix IB] Other extensions, like polygonal trajectories are affected very shortly in Sec. O and 
Appendix [Cl 

2 CONVEX BODY 

2. 1 Basic geometrical models 

There are many different functions and relations between them used for description properties of 
convex bodies [T1[11I31I11[51[S1[71ISIIB EHl Ull 1151 US] ■ In present paper are considered three different 
kinds of distributions: distances between points Fig. [2^, lengths of radii (segments) Fig. [^b, and 
lengths of chords Fig. [2]:. It may be useful to describe precisely models of generation of each 
distribution, to avoid some problems with ambiguity, similar with widely known Bertrand paradox 

mm- 

Definition 1. Distances distribution function is Fri{l) — ri{x)dx with density r] (I). The distances 
are defined by pairs of points inside of the body 5J. Both points are from independent uniform 
distributions, Fig. 

Definition 2. Radii distribution function is Fi^{l) — i{x)dx with density The radii are 

defined as segments of rays from a point inside of the body 23 to the surface. The points are from 
uniform distribution and directions of the rays are isotropic, Fig. O). 

Definition 3. Chord lengths distribution function is Ff^^l) = ^{x)dx with density The 
chords are defined by intersection of the body 2J with isotropic uniform distribution of lines, Fig.\3p. 




Figure 2: Distributions: a) Distances between points, b) Radii, c) Chords. 

It is also convenient to use autocorrelation function j{l) related with distances distribution for 
three-dimenional body as 

rj{l) = Anl'^-fil) /V, dMl in Appendix © 

where V is volume of 03. The expression for autocorrelation function for body with arbitrary density 
Eq. (lASP together with derivation of Eq. (jASP for constant density is recollected for completeness 
below in Appendix IA-21 

There is remarkable correspondence between these densities [TJ [H [31 131 [SI [U [TSl [IS] : 

i^MO = -^'(0 = 7"(0, (2.1) 
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where average chord length (/) ~ lfi{l)dl may be expressed via volume V and surface area S 
using a relation for three-dimensional convex body derived in XIX century by Cauchy, Czuber and 
rediscovered later by Dirac et al ^ [TUl [HI [H] 

(0=4|. (2.2) 

Distribution of lengths of radii [16 is also known as interior source randomness [15j . Relation 
between /i(Z) and in Eq. (|2.1|) is often represented in integral form [151 fTB] 



fi{x)dx^{l)-\l~F^{l)). 

Proportionality between /i(Z) and second derivative of autocorrelation function 7(/) in Eq. (|2.ip is 
also well known and widely used [TJ [H [51 [5] . 

Here is convenient for completeness of presentation and further explanation of nonconvex case 
to derive all equalities in Eq. ()2.ip . considering and r]{l) as generalized functions. 

The usual definition of generalized function [12] is a continuous linear functional T{(j)) on a 
space of a test functions <j). For an integrable function 7p the functional is defined for a test 
function 4>{x) as 



T^{^) - / i;{x)4>{x)dx. (2.3) 

^ — OO 

The generalized derivative [12j is defined as functional 

T'{4>) = ~T{<I>'). (2.4) 

Such a definition ensures derivative of any order for with arbitrary integrable function ^ and 
it justifies use of generalized functions and derivatives in Eq. (|2.ip . 

2.2 Dirac's method of chords 

The Dirac's method of chords |14] uses transition from six-dimensional integral over pair of points 
in some convex body QJ to expressions with chord length distribution, e.g.^ 

V^{^) f f ^drdr' = A /°°^(^) ( rr^{r)drdp] dx, (2.5) 



V" 7<H Jot 47ri?2 4]/ 

where R = |r' — r| is distance between points, dr and dr' are two three-dimensional volume 
elements, S is surface area, and V is volume of 9J. Such an equation was derived in Ref. [14j for 
particular function ip{R) — exp(— ai?). 

It is shown in Appendix [^ that derivation of Eq. (j2.5p has direct connection with equalities 
Eq. p.ip . Here is only outlined basic results and few steps of derivation Eq. (|2.5p . 



First, the integral may be expressed via function of distances ri{l). From Eq. (|A2|1 with 
Eq. (jA4p follows quite understanding relation 

=111 ^g**' = r^.,i^)<i^- (2-6) 



Due to Eq. (|A7|) with Eq. (|A4p right-hand side of Eq. (|2.6p may be rewritten using autocorre- 
lation function 7(;) Eq. (|A5P 

/■OC 

25aj(¥') = / 7(a;)V'(a;)(^a; (2-7) 



© A.Yu.Vtaso'U, 2007 



4 



Signed CLD 



The six-dimensional integral P^j may be reduced to a four-dimensional one by two steps [14] 
revisited in details in Appendix [X] In accordance with Lemma [I] in Appendix lA-ll Eq. (|A9|) and 
Eq. (|A4p it may be expressed 



'Dmi.v) = ^JjrJ dnj^ "'^^ip{R)dR, (2.8) 

where Rmax is length of radius for given point and direction, dr dQ is five-dimensional integration 
on all points and directions. Next, Lemma [H in Appendix IA-31 let us rewrite Eq. (|2.8p using radii 
density function 



V<a{ip) = J l{x) ip{r)drj dx. (2.9) 

Finally, due to rather standard arguments IT^ , revisited in Appendix IA-5[ it is possible to 
rewrite Eq. (|2.8[) using four-dimensional integration on space of lines T 



dp / 

"'0 



2?2t(^) = ^ l^-r I dp I ^{r)dr, (2.10) 



where dT is canonical invariant measure on the space of lines and Lch is length of chord for given 
line intersecting body QJ. 

Now it is possible to rewrite Eq. (|2.10p using Lemma from Appendix IA-51 and Eq. (|A33p to 
produce initial Eq. (|2.5p . 



These integrals also justify use generalized functions, because may be associated with lin- 
ear functionals on some test function ip. The relations between integrals may be considered as 
transformations of these functionals without necessity to indicate any particular (p and it is quite 
reasonable, because main purpose of this paper is rather discussion about geometrical distributions 
than about calculation of some integrals. 

The Eq. (|2.7p may be rewritten using Eq. (|2.3p for space of test functions defined on some 
interval < a; < /max 

V<a = T^ = j. (2.11a) 



The Eq. (j2.9p may be rewritten due to definition of generalized derivative Eq. 

X>2j(0') = / i{x)(l){x)dx ^ V'<^ = -T, = -i. (2.11b) 







Finally, from Eq. (j2.5p follows 

S S SI 

'^■v{(l>") ^ ^ n{x)(t){x)dx =^ V'^ = —T^ = —fi = —fi. (2.11c) 

The Eq. (|2.1ip correspond to Eq. (|2.ip for generalized functions and derivatives. 

The derivation of integral Eq. (|2.7p with autocorrelation function "f{l) is quite straightforward 
and revisited in Appendix IA-21 Eq. (jA7P . The Dirac expression Eq. (|2.5p is directly derived from 
Eq. (|2.7p via two integration by parts if Eq. (|2.ip is true. So if we could consider Eq. (j2.ip as a 
"definition" of a function /i(Z) via second derivatives of 7(/) it ensures Eq. (|2.5p . 

Such a property was used in some works for definition of formal (generalized) chord length 
distribution via Eq. (j2.ip for body with arbitrary shape and density [S] [6] . 
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3 NONCONVEX BODY 



3. 1 Formal integration by parts 

Let us consider application of Eq. (j2.5p to some nonconvex body Ul. The distances distribution 
and autocorrelation function is defined for nonconvex body and so it is possible to write 

oo 

j{x)ip{x)dx (3.1b) 
7'(x) ip{r)dr^ dx (3.1c) 

oc / j^X pp 







7"(x) / Lp{r)drdp]dx. (3. Id) 

\"'0 "'0 / 

Here Eq. p. lap and Eq. (j3.1bp coincide with Eq. (|2.6p and Eq. (|2.7p in convex case respectively, but 
Eq. (|3.1cP and Eq. (j3.1dp are produced by formal integrations by parts and need for geometrical 
interpretation. Let us denote for nonconvex body the signed radii and chord densities represented 
via (normalized) first and second derivative as i±{l) and fJ-±{l) respectively. 

3.2 Radii (signed) density function 

An analogue of Eq. (|2.8p used for introduction of radii density function has form 



im J JRnm 

where R H *yi is intersection of a body 91 with a ray from a point inside the body Fig. [3] 




Figure 3: Scheme of intervals for radii in nonconvex body 
It is possible to write 

/ if{R)dR = XI / v{R)dR, Ro = 0, (3.3) 
JRn<n j=Q-'R2, 

where ni is amount of intervals [R2j, R2j+i\ of ray inside body 91 with R^ for fc > corresponding 
to 2n\ — 1 points of intersection of ray with surface of body Fig. [31 The Eq. p.3p may be formally 
rewritten 

X / ^{R)dR = (-1)'"' / viR)dR. (3.4) 
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Let us denote rimax maximal possible amount of intervals for given body and define ik{l), 
k = 1, . . . , 2nniax — 1 as density function for k-th interval length Rk- Then it is possible to derive 
from Eq. (^21) for 



^±(0= E (3.5) 
fc=i 

analogue of Eq. (12. 9p 



2 r im ax 1 /.QC / /"a 



ip{r)dr j dx — j i±{x) ( / ip{r)dr ) da; (3-6) 

fc=i ^-^0 
It is convenient to consider three distributions '-+(/), and where 

E ^2fe(o, E ^2fe+i(o, i±(o = ^i(o + H(o-^-(o- (3.7) 

fc=i fc=i 

Here "positive" and "negative" distributions, i.e., + and respectively may be 

nonzero for the same Z, e.g., some radii R2 may be equal to R3 in other points. So Eq. ()3.7p could 
not be considered as true Jordan decomposition defined as difference of two nonnegative functions 
with nonoverlapping support [12j 

t±(0 - t"(0, V/ : L+{1) > 0, L-{1) > 0, L+{1) L-{1) ^ 0. (3.8) 

In fact, for some nonconvex bodies l±{1) > and so = despite of nonzero because of 
iiil) + i+{l)>i-il). 

Anyway, the scheme discussed above and expressed by Fig. [31 Eq. (|3.5p . Eq. p.6p . and Eq. p.7p 
provides a stochastic interpretation of i±{l). Similar with Definition\^ there is an uniform distri- 
bution of points inside a body and isotropic rays, considered as some kind of "primary events" for 
radius i?i and Each ray intersecting body *Xl more than one time nj > 1 also produces two 

kinds of "secondary" events: ni — 1 radii i?2fe+i from a "positive" distribution and rij — 1 

radii i?2/c from a "negative" one 

Such a stochastic model also produces understanding description of some integrals via averages, 
mathematical expectations etc., e.g., 



/ Li{l)dl^ / L±{l)dl^l, / i+{l)dl^ / 

JO Jo Jo Jo 



i-{l)dl, (3.9) 



because each radius corresponds to ni "positive" events and ni — 1 "negative" events, i.e., contri- 
bution to "total charge" or "balance" of events N = N+ — iV_ is always one. 

It is also useful to introduce distribution of total length of all segments for given radii to(0 
(OSD, one-segment distribution), then 

{l),= / L±{l)ldl^ / 10(0'^^^ (3.10) 



because contribution for any ray is J2h~i ^^fe = -Ri + X]fcLi^(^2fe+i — R2k), i-e., total 

length of all segments. 
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3.3 Chord length (signed) density function 



An analogue of integral Eq. (|2.10p for nonconvex body 91 may not use continuous area of integration 
on dpdr, if chord intersects 01 along few intervals. It is similar with integration along set of interval 
for radii Eq. (|3.2p and Eq. (|3.3[) in Sec. I3.2[ but decomposition is more complex because of two 
integrals. 

For a chord with rij intervals there are 2ni points of intersection (Lq = 0,Li, . . ., L2k, L2k+ii 
. . .) for fc = 0, . . . , ni — 1. In Eq. (|2.10p p is coordinate along the chord and r is distance between 
points. If to introduce new variables x = x' — p + r, then for whole chord inside body 

A^{(p)= / (p{r)drdp = - / (p{x' - x)dx dx' ^ - / (p{x' - x)dx' dx, (3.11) 
Jo Jo Jo Jo Jo Jx 

but for few segments the integration of ip{r) ~ ip{x' — x) on dxdx' should include all points 
x,x' £ Ln^, X < x', where L n Dl = [jh=o i^^k, L2k+i] is union of all intervals of given chord 
inside the body, Fig. [4^. Here the minus signs in Eq. (|3.11[) are due to term —x in Lp{x' — x). 



^5 X 




u 








L2 










c) 



-A. 



Figure 4: a) Scheme of integration, b) Scheme of intervals c) Decomposition of Bj^fc 

Let us denote (L n Ol)^^ area of integration described above and depicted on Fig. as set of 
triangles and rectangles below of the diagonal x — x' 



rii — 1 k<ni,j<k 

(LnoT)^ = U^'^U U 

fc=0 k=l,j=0 

Afc = {{x, x') : L2k <x<x'< L2k+i}, 

= {{x,x') : L2j < X < L2J+1, L2k < x' < L2k+i}- 



Now it is possible to write analogue of Eq. (|2.10p for nonconvex case 



(3.12a) 

(3.12b) 
(3.12c) 



(3.13) 
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Due to Eq. (|3.12p it may be rewritten 
/ / '^{x' — x)dxdx' = 

= -(^^ 1 1 (p{x' -x)dxdx' 1 1 if(x'-x)dxdx'^, (3.14a) 
where 

tf{x'~x)dxdx'= ip{r)drdp:^ k^'"'+^~^'"'{ip) (3.14b) 



'(Lnm) 

ni — 1 „ „ rii — l k—1 



and integral on Mj^k may be decomposed on four terms illustrated on Fig. [4j: 
ip{x' — x)dx dx' = — I I ip{x' — x)dx' dx 

lj,fc L^j L2k 



2j 



I / Lp{x' — x)dx' — / ip{x' — x)dx' jdx 



L2J + 1 /•L2k rL-2j rL2k (3.14c) 

tf{x' — x)dx'dx — / / Lp{x' — x)dx'dx 



'0 Jx Jo Jx 

<'L2j + l pL2k + l f^2j pL2k+\ 

I ip{x' — x)dx' dx + / (p{x' — x)dx'dx 

Jx Jo Jx 

If ni is number of segments, there are ni ("positive") terms in Eq. (j3.14bp and /^Hiillizll terms 
(nj — ni "positive" and — ni "negative") in Eq. p.l4cp . 
It is possible to decompose Eq. p.l3p using Eq. p.l4p 

1 /"/"'"^ 



nj — 1 k—1 



k=l j=0 

ni — l k~l 

- + ^^"'^'^^'''+'{f)])d'T'. (3.15) 

fe=l j=0 

Let us recall Eq. (|2.10p for convex body expressed using Eq. (j3.14bp 

-D^iv) = ^ / A'-(^)dr. (3.16) 

The Eq. ()3.15p with three terms corresponds to sum of 2nf — ni integrals Eq. p.l6|) arranged in 
three groups. The first one includes all ni segments [L2k, -^2fe+i] of given chord. Second one takes 
into account ni{ni — l)/2 pairs of intervals [L2j, L2k+i] and [L2j+i, L2k]- Third "negative" term is 
for ni{ni - l)/2 pairs of intervals [L2j,L2k] and [L2j+i, ^2/0+1]- 

For a chord with ni intervals there are nj "positive" and — ni "negative" terms. The scheme 
is depicted on Fig. [IJj for m = 3, there "negaitive" intervals are drawn by dashed lines. 
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It is possible to decompose ^i±{l) on three parts corresponding terms in Eq. (|3.16p : 

fj.+ {l), fJ'-il)- Here the fii{l) up to normahzing muhipher corresponds to density /im of muhi-chord 
distribution (MCD) discussed in Sec. [TJ 

/■oo 

Aii(0 = CMm(0, CM = / ni{l)dl. (3.17) 
If normahzation of fJ-±{l) is required, it should be expressed 

M±(0 = Cm'[mi(0 + M+(0 - A*-(0] = mil) + cm'[m+(0 - M-(0]- (3.18) 

Here again "positive" c^[fii{l) + fJ,+ {l)] and "negative" c^^i-{l) terms of Eq. p.lSp can be 
overlapped and so formally could not be considered as Jordan decomposition [12j 

Ai±(0 = -"+(0 - -""(0, V/ ; > 0, > 0, = 0. (3.19) 



Let us write finally analogue of Eq. (|2.5p for nonconvex case 

V<n{^) ^^J^J^ dr' - i^' ^"^^±(2:) (^Jjyir)drdp^ dx, (3.20) 

where ityi is some constant. It was shown earlier for convex body — (l) = 4:V/S. 

Definition of signed chord distribution as quantity proportional to second derivative of auto- 
correlation function suggested in Sec. [T] together with comparison of Eq. 13:20)) and Eq. ((3ld|) 
produces rather formal equations like 

^mV±(0 = 7"(0 = -4(0, - r ^"(Ocf^ = 7'(0), (3.21a) 

but Eq. (|3.18p also describes "formal" /i±(0 via geometrical distributions of segments of chord in 
nonconvex body. 

It is also possible to write analogue Eq. (jA32[) 

/•oo /'OO poo 

£^=i<j^ L±{l)dl = -e<jx L'^{l)ldl^ I fl±{l)ldl ^ (l). (3.21b) 
Jo Jo Jo 

If to use Eq. (|3.20p with (p{l) = AttP it is possible to derive from Eq. (|3.1d[) yet another one relation 
for an average 



dr dr' _ V"^ _ f°° ^ fj^^TTl"^^, , _ 7r(/'* 



y -y=^.'l M±(0^rf^ =^ '- = iv^ (3.21c) 



It is possible to use stochastic model to clarify some expressions above and derive other useful 
results. There is uniform isotropic distribution of lines described in Defimtion\^aiid corresponding 
to "primary" events. If such line intersects body nj times, there are 2nf — ni "secondary" events. 
There are n/ segments of given chord together with ni{ni — 1) "positive" and "negative" intervals 
described above. 



Note on event counting: An essential difference with radii distribution is contribution ni > 1 
for each chord to "total charge" iV = — N^. So number of lines Ni is not equal with TV and 
Ni/N Cm for N — > 00. The similar effect is true for work with usual MCD case and so should 
not be considered as specific difficulty of signed chord distribution. 
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It is also possible to express cm using /io(0 density for one-chord case (OCD) also mentioned 
in Sec. [1] because contribution to average length for fii (without normalization on cm) is sum of 
all segments and it is the same for OCD case 

xD noo /"Oo /"Oo n\ 

lio{l)ldl= / ^i{l)ldl — CM ^J,Mil)l dl — Cm fi±{l)ldl, Cm = -ttt-- (3.22) 



(0 

There is also quite similar expression for (P). The contribution to P for one chord due to 

Eq. 

ni — 1 

^2fc)^ 



{L2k+1 — L2 
11-1 k-1 

^^((i2fc — i2j + l)^ + {L2k+1 — L2j)'^ — {L2k — L2j)'^ — (L2fe+1 — L2j+lY 



k=0 

m — l k—l 



2 



fc=l j=Q 

ni — 1 ni — 1 k—l 

= ^ {L2k+1 — L2k)^ + 2 ^ ^(i2fc+l — L2k){L2j+l — L2j) 
k=0 k=l j=0 

ni — 1 

= (5I(Wi--^2fc))'' (3.23) 

fc=0 

Due to Eq. p.23p this contribution (to fii + /i+ — M-) is always equal to square of sum of all 
segments, i.e., coincides with OCD case. So, due to Eq. (|3.23p and Eq. (|3.18p 

XD POO POO //2\ 

tio{l)l^dl^ [fii{l)+fi+{l)-fi-{l)]fdl^cM fi±{l)l^dl, CM = ^. (3.24) 



Finally, let us recall few methods for calculation of multiplier £<yi. It is discussed in Ref. [7] 
that for convex and nonconvex case 7'(0) — S/AV and so due to Eq. (I3.21a|) . it is possible to 
use £<ji — 7'(0)^^ — 4:V/S similar with nonconvex case. There are also interesting results [U [2] 
concerning {l)o and {1)m for nonconvex body. It was proven for wide class of nonconvex bodies [5] 
{1)m = 4V/S, and so (l) = {1)m = -iV/S and Eq. ((XITbll again ensures e<yi = (l) = AV/S. 

It was also mentioned [2] yet another useful result (/)o — AV/S*, where S* is surface area 
of convex hull and so due to Eq. (|3.22p cm = S/S*. On the other hand, it is not quite clear, if 
equations with S* are true for nonconvex body that may not be represented as convex body with 
few convex holes. 

The stochastic model considered here also clarify use Monte-Carlo methods for calculation of 
signed chord length distribution used in Eq. (|3.20p . It is necessary to generate uniform isotropic 
random lines using methods discussed elsewhere [T7]. For each chord, intersecting a body ni times, 
lengths of all 2nf — n\ segments represented in Eq. p.lSp should be taken into account with proper 
signs. 

It is possible to generate either one signed distribution or two distributions for "positive" and 
"negative" segments to represent n±{l) via difference like Eq. (|3.18|) . It is also necessary to use 
proper normalization, i.e., use the "total charge" N = — 7V_, rather than number of lines Ni 
(see Note on event counting above on pagellOp. 

4 NONUNIFORM CASE 

The nonuniform case has more ambiguity in definitions and difficulty with geometrical interpreta- 
tion. It is only briefly mentioned here. Two analogues of Eq. (|2.5p are considered in Appendix |B] 
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and have form 



J I' pir)p{r')^^^drdr' = C, jy^{x) (^Jjj{r)drdp^ dx, (4.1) 

where Ary is some function, Ca is constant and ij-a{x) is a density for chosen function Ar.r'- 

A simple choice Ar.r' = \r — r'\ is discussed in Appendix IB- II Eg. (jBip and it is proven the 
correspondence with definition of CLD ^ /^a(0 via second derivative 0[6]. 

In some apphcations it is more appropriate to consider "optical" length Ar^ — O'^ defined by 
Eq. (jB4p as an amount of substance along straight line between r and r' . For uniform case and 
convex body it is proportional to distance between points, but for nonuniform case it is another 
integral Eq. (|B5|) reproduced in Appendix IB-21 

For such an integral there is also analogue of Eq. (|2.5p . It is expression Eq. (|B13|) with distri- 
bution of "optical" chord lengths fl{l) ^ IJ^t\(J) used in Lemma in Appendix IB-21 The complete 
proof of Eq. (jB13[) with formal "cancellation" of two terms with p may be found in Appendix IB-21 
This case may be considered as even more direct analogue of convex body with uniform density 
Eq. p.Sp . because is certainly nonnegative due to quite clear geometrical definition. 

The p{l) is also may be used for nonconvex body *n with uniform density p = 1 by formal 
consideration of union of given body and complement to convex hull with p — Q. In such a case 
"optical" length of chord corresponds to sum of all intervals inside 01, i.e., definition of OCD case 
[H m [3] mentioned in Sec. [1] 

The application of OCD case for nonconvex body with uniform density was suggested already 
in initial paper about method of chords 14J. 

5 APPLICATIONS TO ARBITRARY PATHS 

A natural extension of problems discussed in present paper — is the consideration of polygonal 
paths. There are interesting results, concerning trajectories of random walk, e.g., validity of Cauchy 
formula Eq. (12. 2p for average path length inside the body ^ [TH H EH |22] . Appendix [C] contains 
a rather illustrative geometrical derivation of this result for uniform isotropic case. 

Possibility of application to some tasks concerning arbitrary trajectories is important property 
of method of chords both for convex Eq. (|2.5p and nonconvex case Eq. (I3.20p . Sometimes consid- 
eration of such integrals uses understanding, but simplified models with propagation of particles 
along straight lines. Such a suggestion is not always justified, because in most cases considered 
above the segments, radii and chords — are formal subspaces of integration. The only necessary 
condition is isotropy and uniformity of model, justifying use of function (p(|r' — r |), instead of some 
general function if{r' ,r). 

Formally, the only example of more general function mentioned here is ip(0^?) used in Eq. (|B13|) 
and discussed in Appendix IB-21 For such a model consideration with propagation along line is 
justified due to definition of O^, in Sec. Hand Appendix IB^ Eq. (|B4|. 

Here should be compared two cases. One model with propagation along lines makes possible 
to use method of chord for convex, nonconvex and even nonuniform case Sec. |4] and Appendix 
IB-21 Other model may be used only for uniform isotropic case, but it is possible to consider also 
propagation with scattering, that may be described by some function Lp{\r' — r\) due to spherical 
and translational symmetry Fig. [5l 

The basic theme of presented paper is the second model described by signed density function 
H±(l). The first model was discussed for completeness in Sec. |4] and Appendix IB-21 and it is shown, 
that it always described by some nonnegative density function without essential difference in 
construction for nonconvex and even nonuniform case. The construction of signed chord density 
function is much complicated and was described in Sec. 13.31 
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Figure 5: Trajectories and lines: a) Nonconvex body b) Convex body c) Spherical symmetry in 
distributions of trajectories 

Let us consider examples of application of both cases. On Fig. [5] is represented propagation 
along some trajectories. For nonconvex case Fig. [5^ any trajectory may intersect body more than 
once. For straight line it is always possible to use first model with For nonconvex case and 

environment with neglectful density OJi, is simply sum of all segments inside body along given line. 

A polygonal trajectory may intersect body more than once for both nonconvex and convex 
cases Fig. [S] If body and environment is the same uniform isotropic media, it is possible to use 
integrals with </j(|r' — r|), otherwise there are "boundary effects," because different trajectories 
intersect environment by different ways depending on concrete positions r' and r even for convex 
case Fig. [5)d. 

An illustrative example of an application is a set of random trajectories from uniformly dis- 
tributed points r inside the body with isotropic initial directions (see Fig. [5t) and function 
= y(|r' — r|), treated as probability density of cancellation of trajectory (absorption) on 
distance R from an origin r. Then Eq. (12. 5p or Eq. ()3.20p describe probability for a trajectory of 
such random walk to be finished inside the body. 
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APPENDIX A CALCULATION OF DISTRIBUTIONS FOR CONVEX BODY 



A-1 Distribution of distances 

Let us consider some convex body 2J C M'^, function $(a;) and integral 

I<a{<^) ^ I I <P{\r' -r\)drdr', r^{x,y,z), dr^dxdydz, (Al) 

where V is volume of 03. It is six-dimensional integral with the function $ of a distance |r' — r| for 
pair of points r, r' G 5J. The multiplier is used for normalization X<jj(l) — 1 and Eq. (jAip may 
be also treated as an average of the function $ of one variable and expressed via single integration 
due to lemma below. 

Lemma 1. The linear functionaU<s{^) defined by Eq. \Al\l may be rewritten 

poo 

X<h($) = / Hx)rj{x)dx, (A2) 







where ri{x) is density of distances distribution in body 5J introduced in Definition]^ 

Proof. The linear functional I(jj(<i>) also should be considered as a generalized function and for 
regular case it may be represented via integral like Eq. (|2.3p . i.e., Eq. (jA2p with some r]{x). Let 
us show, that r]{x) is density for distribution of distances. 

It is enough to consider Iqj with functions ^{x) — Q{1 — x) = Qj'{x), there 8 is Heavyside 
step function, i.e., Qj'^x) is equal to zero for x > I and unit otherwise. For such a function Xgj 
integrates over all pairs of points in 03 x 9J with distances less than I and due to Eq. (|A2ll it is 
possible to write equation for probability X(h(8^) = P(|'''' ~ < l) — Jq ri{x)dx = F,,(0 and it 
coincides with definition of density ri[x) for distances distribution function Fri{x). 

For extension of the proof for non-regular case, it is possible to use directly the distances 
distribution function F,i{x) and to write instead of Eq. ()A2|) Lebesgue-Stieltjes integral [121 123] 

/•oo 

= / $(x)dF^(x). (A3) 

□ 







For consideration of both regular and non-regular cases it is also possible to treat as a 
generalized function. In such a case representation of T<s{^) as an integral Eq. (jA2p may look as 
a not very rigor representation of a tautology like Xfjj = 77 = T,, , because a generalized function by 
definition is a linear functional [12] . 

The functional Xqj is connected with Pqj defined in Sec. 12.21 Eq. (|2.5p by straightforward 
relations 

Vv{^)^VI<n{^), I^{<S>)^^V^{47rxH{x)), (A4) 

in( X) 

produced by the choise <i>(a;) = V- — tt. 

Attx'^ 

The meaningfulncss Lemma[I]is not only integral representation Eq. (jA2p . but association of Xqj 
with distribution of distances. It is simple demonstration of some methods to avoid Bertrand-like 
paradoxes, because ri{x) is not only possible density for distribution of distances between points in 
a body. 
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For calculation of Ifjj Eq. (|A1[) is used integration on six-dimensional space with natural Eu- 
clidean measure on = R'^ x M.^ and it is in agreement with independent uniform distributions 
of both points used in Definition [7] 

An example of alternative measure was used in Ref. [THj: it was considered distribution of 
chords like in Definition above and segment of line between pair of points on the chords. A 
distribution function for lengths of the segments, i.e., distances between ending points may be also 
calculated using Heavyside step function [16j . but density is not equivalent with ri{x). 

Yet another example of similar measure may be constructed if we consider first point r from 
uniform distribution, but the second one r' is generated with isotropic distribution of relative 
directions R=r' — r and uniform distribution of distances \R\ = |r' — r\ between points. In such 
a case distributions of points r and r' are correlated and relatively to natural Euclidean measure 
on {r,r') g R'^ x R'^ = R^ D 9J x 5J here is necessary to introduce multiplier |-R|^^ — \r' — r\~^. 

Up to constant normalizing multiplier the alternative density of distances introduced above 
may be expressed as r]{x)/x'^, i.e., proportional to autocorrelation function [U 31 [31 [B] discussed 
further. 

A-2 Autocorrelation function 

For body with density p{r), r G R'^ is defined the autocorrelation function 7(r), r € R'^ or 7(/), 



i.e., 7(Z) is an average of 7(r) on sphere with radius /, {Si : \r\ = I}. 

In simplest case of constant density p{r) = 1 for r G 2J and zero otherwise. For such a case 
7(0) = V, and it is convenient to rescale g{r) = p{r)/y/V for normalization 7(0) = 1. It is possible 
to rewrite Eq. (|Aip 



where Eq. (jA7p is produced from Eq. (jA6p by integration over spheres S;. 

Comparison of Eq. (jA7p and Eq. (jA2p with arbitrary function produces relation between 



I G M 




(A5) 




(A6) 



(A7) 



and 



V{1) = TTiHl) 



(A8) 



already mentioned earlier in Sec. 12.11 
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A-3 Distribution of radii 



Let us introduce new variable R — r' ~ r in Eq. (jAip and rewrite integral on dR using spherical 
coordinates 

1 r /'27r pR(r^9,4>) 

= — dr sinOdO / d0 / R^^{R)dR 
^ J<n Jo JQ JQ 

= —^ dr dfl R^^{R)dR, (A9) 

J<n Js Jo 

where R^ sinO is Jacobian in spherical coordinates (R, 9, (p) of vector r' — r, § is unit sphere, and 
R{r,9,4)) — R(r,r2) is length of radius (segment) mentioned in Definition\^ i.e., distance from 
point r to surface of body in direction represented by spherical angles 9, (j) or unit vector G §, 
Fig. El 




Figure 6: Scheme of integration along radii with length R(r, 9, (p) = \R\ 
Let us introduce linear functional 

where r e C M'', e §, and R(r,r2) is notation for length of radius introduced in Eq. (jA9p . 
Here ~ 1 and AttV is normalization, i.e., 5D volume of set QJ x § with respect to canonical 

measure on R"^ x S. It is analogue of normalization of Eq. (|A1[) with V'^, i.e., 6D volume of set 
2J X 5J with respect to canonical measure on R'^ x K^. 

Lemma 2. The linear functional defined by Eq. iA 1 0|) may be rewritten 

^^\'^) = / '^ix)i{x)dx (All) 
Jo 

where l{x) is density of radii distribution in body QJ introduced in Definitions^ 

Proof. The proof is similar with Lemma{^ The Eq. (IA11|) is again almost tautology for generalized 
functions zj^'' = l = T^. Let us anew consider I^"* with Heavyside step function 0^(a;), to integrate 
over points in 5J x § associated with radii less than given /. 

The measures of integration in correspond to uniform distribution of points r and isotropic 
distribution of directions Jl on unit spheres (due to transition to spherical coordinates in second 
integral). It is in agreement with Definition\^and soI^^(6^) = P(R(r, fl) < I) — i{x)dx = F^{1) 
and meets definition of density t(x) for radii length distribution function i^i(x). □ 
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The functional J^^ Eq. (|JTO| let us simplify Eq. (|M 



47r r 



= (^y R'<P{R)dRj (A12) 
and substition of Eq. (|A12p to Eq. (jAlip produces yet another expression for Xqj (<i>) 



X^{<^)^ j i—j R^^{R)dRy{x)dx. (A13) 

A-4 Autocorrelation function and distribution of radii 

Due to definition Eq. (|A8p it is possible to rewrite Eq. (jA2p 

poo A 

Jaj($) = / ^{x)—x^-i{x)dx. (A14) 

Yet another expression 

(■OO / A pX \ 

Imm = -J iy J R^^{R)dRU'{x)dx, (A15) 

may be produced from Eq. (|A14[) using integration by parts. The equality of Eq. (jAlsp and 
Eq. (|A13P corresponds to second equation in Eq. (|2.1|) 



- l{x) = i{x). (A16) 

The integral representations above let us also consider Eq. (jA16P as derivative of generalized 
function. 

A-5 Ctiord iengtti distribution 

Let us change order of integration in Eq. (|A9|) 

^ dfl dr (A17) 

V Js Jo 

and for any unit vector O e § to decompose spatial integral J dr — J J dP dn on two integrals: 
along the axis n — no parallel to Vl and on the plane P = Pq perpendicular to ft Fig. [71 




Figure 7: Lines and chords 
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After such decomposition Eq. (jA17|) may be rewritten as 

If r /•L(ps!,n) 
^ dn dP dn R'^<i>{R)dR, (A18) 

^ Js JPnCX}) Jo Jo 

where Pn{^) is projection of 03 on plane Pn and L(pn,r2) is length of chord of line defined by 
direction fJ S § and point pn G Pn, see Fig. [T] 

It is possible to define four-dimensional space T of (directed) lines [51 [21 [13 [HI [13 US] , there 
each line is defined by point on sphere S S and point of intersection pQ with orthogonal plane 
Pq _L n. Here each line is represented twice with two opposite directions ±d. It is possible to 
introduce space of undirected lines as a quotient space T = T/j+l, —1}, but in most expressions 
below for convenience of calculations is used T. The first two integrals in Eq. (jAlSP correspond to 
integration on this space 



l%^{<^) ^tt:^ dT dL R^^{R)dR, 

^ JT{<a) Jo Jo 



(A19) 



where T(53) is set of lines intersecting 53, \-{It) is chord length for a line Ir G T , and dT is canonical 
(uniform and isotropic) measure on T unique defined up to constant multiplier by invariance with 
respect to translations and rotations [H [3 [T3 [26] . 
Let us introduce linear functional 

4'^(T) = TTy^, I T[L(/r)]dr, V[Tm = / dT, (A20) 

where y[T(53)] is normalization, i.e., 4D volume of T(23) with respect to dT . It may be written 
also due to definition of J.j-^^^dT above 

V[Tm = jdnj dP^ ^<Sp,,(^)) (A21) 

where {Sp^(<}}-yi is average surface of projection of body 53. For convex bodies due to a Cauchy 
formula [311] 

(5p,(«)) = (A22) 



where 5* is surface area of 03 and so 



V[T{^)] = 7r5. (A23) 



Lemma 3. The linear functional defined by Eq. IIA20\) may be rewritten 

poo 

J^''(T)== / T{x)fi{x)dx (A24) 
Jo 

where fJ.{x) is density of chord length distribution in body 53 introduced in Definition\^ 

Proof. Similar with Lemma [7] and Lemma [^ the Eq. (|A24p may be represented in rather trivial 
form I^'"* = ^ = Tf^ for generalized functions. Here again the Heavyside step function 6^(a;) selects 
subset in T(93) associated with chords shorter than given L. The measure of integration dT in 
2'^''' corresponds to uniform and isotropic distribution of lines It in agreement with Definition\3l 

and so X^'^(8^) = P(L(/7-) < L) = fi{x)dx = Ff^{l) and meets definition of density /i(x) for 
chord length distribution fmiction Ffj_(x). □ 
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The space of lines has nontrivial structure and an essential moment in proof of Lemma is 
correspondence of measure of integration dT in Eq. (|A20[) and measure describing distribution of 
lines in Definition This isotropic uniform distribution sometimes associated with concept of 
/Lt-randomness [Tl [T?l [TO] . 

An alternative distribution of lines is often defined by uniform distribution of points inside the 
body 03 and isotropic distribution of directions. It is sometimes called i/-randomness [TJ I19j or 
interior radiator randomness [15j . It is different from /i-randomness, because the same line may 
be represented by any point on this line and direction along this line. So measure of i/-chord for 
given line in comparison with uniform case has extra multiplier proportional to length of the chord 
v{l) oc The precise expression for convex body is pi ITSl fT9] 



KO = ^MO=(^^j (A25) 

The concept of i/-randomness for chords is close with Definition [H of radii distribution and it is 
useful for some applications [TH [12] . 

Yet another distribution is produced by definition of line by pair of points with independent 
uniform distributions inside 03. It is sometimes called A-randomness [1] [6l [19] . The measure for 
A-chord for three-dimensional case A(^) oc The normalizing multiplier for may be directly 

calculated [H H [TOl [H [S] 

= ^ (A26) 
and precise expression for convex body is [TJ (TH] 

A(0 = ^m(0 = i^-^j I'Kl)- (A27) 

The concept of A-randomness for chords has certain relation with Definition [7] of distances distri- 
bution. 

Such different kinds of randomness illustrates necessity of rather pedantic work with distribution 
of lines due to analogues of Bertrand paradox [8l I13j already mentioned earlier. 

Lemma 4. The functional Eq. iA20\) may be formally expressed via Eq. iAlO\) 



l'^^il)dl]^^J^\^). (A28) 



Proof. The derivation of Eq. (|A28p is similar with transition from Eq. (|A9|) to Eq. (jA19|) via 
Eq. ^K^7\ and Eq. (lUSll 



S AttV J<s Js Js J<n 

\ f f /■L(pn,n) pn / \ / 

—^ dfl dP I dn / 

T^S Js JPn(aj) Jo Jo \Jo 



□ 



An application Eq. (|A24ll and Eq. (|A11|1 to Eq. (|A28j) produces 

^ (^J^ ^{l)dl^ fi{x)dx ^{x)t{x)dx. (A29) 



© A.Yu.Vtaso-u, 2007 



19 



Signed CLD 



An integration by parts of Eq. (jA29p produces 

^ (^J^ ^{l)dl^ H{x)dx = (^J^ ■i'iljdl^ L'{x)dx. (A30) 
The first equality in Eq. (P7T|) follows from Eq. (|A30p 



4:V\ \ , 1 



i'W = -(^^j m(^) = -^A^(^)- (A31) 

It may be considered as generalized derivative due to Eq. (|A29p . 

In Eq. (jA31|) was used Cauchy relation Eq. (|2.2p iV/S = (I). It is possible also to find 
normalizing multiplier for n^ii{x) = ~-l'{x) simply using integration by parts 



i{l)dl = -n~^ L\l)ldl^ lfi{l)dl ^ (1). (A32) 



Jo 

This derivation uses only the fact of proportionality of ij.{x) to derivative of another density function 
and in further applications for nonconvex case normalization with (l) may be preferable due to 
nontrivial proof of expression with volume and surface area. 

Lemma 5. The functional Xs:^ Eq. IIA1\) may be expressed vial'^^ Eq. IIA20\) 

S rf f'^TT 



= ^4" \^J^[j^ yrH{r)drj dlj . (A33) 

Proof. The expression Eq. (|A33|) is straightforward combination of Eq. (jA12p and Eq. (jA28p . □ 
Application of Eq. (IA24P to Eq. (jA33P produces 

The Eq. (Ell) and Eq. (IXM produce for $(r) = V"^^^ 
in agreement with Eq. I|2.5p in Sec. 12.21 It is also convenient sometimes to use expression 



JO 



ip{r)dr dp — / {x — r)ip{r)dr. (A36) 



APPENDIX B SOME EQUATIONS FOR NONUNIFORM CASE 

Let us consider a body with nonuniform density p{r). Any such body ^ may be treated as a 
convex one without lost of generality by consideration of the convex hull QJ and assignment of zero 
density to the complement 93 \ OT. 

Here is considered two cases corresponding to choice of A^-.r' in Eq. (14. ip . 
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B-1 Distance between points Ar,r 
A direct analogue of Eq. (|A7|) is 



, ip(\r'-r\) , f°° 
Pir)p{r ) ^^|^, _^|2 C^^^»^ = j{x)ip{x)dx 

7"(x) i^J J ip{r)drdpj dx, (Bl) 

there second equality produced by two integrations by parts. So up to normalization with 
"f"{x)dx = 7'(0) it is possible to use a formal ("generalized" [3 [6]) chord length distribu- 
tion fi{l) — 7"(/)/7'(0) for calculation of integrals like Eq. (|B1|1 . 
Finally 

j I p(r)p(r')|^^!^^rfr dr' = J^ix) (^ljj{r)drdp^ dx, = 7'(0). (B2) 
There is yet another way to express C^. If to consider (p{l) — AttP, then from Eq. (|B2p follows 



= JJ p{r)p{r')drdr' ^C^ i1{x) (^J^J^nr^drdp^ dx ^ ^C^ 



x*il{x)dx. 



where M = J p{r)dr is mass of the body. Let us also denote (Z^) = x'^il{x)dx, then 

= (B3) 

For case of constant unit density M = V and due to Eq. (|M6| 12^^2/(^5-) (j^ ^ 5/4 

in agreement with Eq. (|A35|) and Cauchy formula Eq. (|A22|) . 

B-2 "Opticai" iengtti Ar,r' = 

In some applications instead of distance between points R= \r' — r\ in tf{R) it is necessary to use 
an "optical width", i.e., integral on density along line between points 

O^' ^}pd£^\r' -r\ [ p{r + {r' - r)x)dx (B4) 

r Jo 



J{^) ^ II p{r)p{r') J^l "U ^rdr' (B5) 



and to consider functional 

47r|r' — r\ 

Let us introduce new variable and rewrite integral on dR using spherical coordinates like in Eq. (|A9[) 

J{ip) ^J-[dr[ dQpir) [ ^ \{r + xn)ip{0^+'=")dx. (B6) 

where R(r, $7) is lenght of (maximal) radius, Ol^^^^ = p{r + lO,)dl = s{x), s'{x) = p{r + 2:0). 
Last integral has form s' (x) p{s{x))dx and may be rewritten as /q"''^ p{s)ds. 



1 



Li(r,0) 



Ji(p) = ^ I dr I dnp{r) I ip{w)dw, (B7) 
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where w{r, n) = Q^^^^'''"''" ^ lo''''^"^ Pi^ + ^^)'^^ "optical" radius length. 

It is possible to use analogue of integration in Eq. (lAlSp with axis n || $7 and plane P ± Jl 

I r r r rw(r,n) 
J{(p)^— dfl dP dnp{r) (p{w)dw, r^rp + nn, (B8) 

where point rp is intersection (projection) of line representing path of integration with plane P 
(see Fig. [7]) and range of n corresponds to variation of r along whole chord in convex span of a 
body. 

If to consider last integral in Eq. (|B8p as some function F{w{n)), two last integrals have form 
/ w'{n)F{w{n))dn = / F{w)dw and so we have analogue of Eq. (|A19p 



1 



r 



J{ip) ^ — I dT I dw I ip{x)dx, (B9) 



where \N{lr) = O^i^'^^K^r) is "optical" chord length. 

Lemma 6. Let us now introduce "optical" chord length distribution fl{x), i.e., to any chord of 
line intersecting (convex hull of) a body in two points ri and r2 instead of \r2 — J^il is assigned 
"optical" length O^J. The integral Eq. fB5\} may be expressed 



poo / rx np \ 

J{(p) = C^J ^i{x) i^J J (p{r)drdpj dx. (BIO) 

Proof. The Eq. (jBlOp follows from Eq. (jB9p and it is complete analogue of derivation Eq. (jA34p 
from Eq. (|A19p in Appendix 1X1 In both cases is used the same integral of some function F{It) 
over space of lines T with the same measure dT and particular method of calculation of F does 
not matter. The fi{l) in Eq. (|A34p is (joint) density for F(?t) = L(/r) and the //(/) in Eq. (|B10p is 
(joint) density for F{It) = W(Zr)- 

There is unessential difficulty with constant multiplier C^, due to lack of simple analogue 
of Cauchy formula Eq. (jA22[) for average surface used in normalization of the integral on T in 
Eq. (|A23|) . An alternative way of calculation of is represented below. □ 



Let us introduce quantity 



II i^^drdr'^ I ,ix)dx (Bll) 



then for Eq. (|BT0)) with (p{l) = 1 

G = C^J^ Kx) (^JJ^drdp^ dx = ^ x^fi{x)dx = ^Cf,{fi) 



and 



C'm = ^ (B12) 



Finally 



Pir)p(r')^0^drdr' = j^fi{x) (^JJ'^ir)drdp^ dx, = ^. (B13) 



© A.Yu.Vtasou, 2007 



22 



Signed CLD 



APPENDIX C AVERAGE PATH LENGTH 



Here is presented simple geometrical proof of equality of average path length inside a body and 
average chord length in uniform isotropic case. 

On the Fig. [5^ is depicted some path ABC with one kink inside a body. Let us also consider 
path FBD produced from ABC by central symmetry with respect to point B, Fig. [Sja. The sum 
of lengths of two paths FBD and ABC is \ AB\ + \BC\ + \FB\ + \BD\ and coincides with sum for 
two chords AF and CD Fig. [St. 




Figure 8: Illustration for average path length 

This method let us get rid of one kink and after few such steps to consider average length of 
straight lines instead of paths. The necessary condition for such a proof — is isotropic and uniform 
distribution of paths to ensure equal probability (density) of ABC and FBD. 
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